Magnetic order, Bose-Einstein condensation, and superfluidity 
in a bosonic t-J model of CP^ spinous and doped Higgs holons 
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We study the three-dimensional U(l) lattice gauge theory of a CP^ spinon (Schwinger boson) 
field and a Higgs field. It is a bosonic t-J model in slave-particle representation, describing the 
antiferromagnetic (AF) Heisenberg spin model with doped bosonic holes expressed by the Higgs field. 
The spinon coupling term of the action favors AF long-range order, whereas the holon hopping term 
in the ferromagnetic channel favors Bose-Einstein condensation (BEC) of holons. We investigate 
the phase structure by means of Monte-Carlo simulations and study an interplay of AF order and 
BEC of holes. We consider the two variations of the model; (i) the three-dimensional model at finite 
temperatures, and (ii) the two-dimensional model at vanishing temperature. In the model (i) we 
find that the AF order and BEC coexist at low temperatures and certain hole concentrations. In 
the model (ii), by varying the hole concentration and the stiffness of AF spin coupling, we find a 
phase diagram similar to the model (i) . Implications of the results to systems of cold atoms and the 
fermionic t-J model of strongly-correlated electrons are discussed. 

PACS numbers: 37.10.Jk, 74.72.-h, 75.50.Ee, ll.15.Ha 



I. INTRODUCTION 

In recent years, cold atoms have attracted interest 
of many condensed-matter physicists. Systems of cold 
atoms have exhibited (and shall exhibit) various inter- 
esting properties like Bose-Einstein condensation (BEC), 
superfluidity (SF) caused by the BEC, magnetic ordering, 
etc. Both for condensed-matter experimentalists and the- 
orists cold atoms offer an ideal testing ground to develop 
and check their ideas because one can precisely control 
parameters characterizing these systems like dimension- 
ality of the system, strength of interaction among atoms, 
concentration of atoms, etc. An example of these ideas 
is a possible interplay of magnetic ordering and BEC (or 
equivalently SF), although coexistence of these two or- 
ders seems not to have been reported experimentally so 
far. 

A standard model of cold bosonic atoms with repul- 
sive interactions may be the bosonic Hubbard model in 
which electrons are replaced by bosonic atoms. For defi- 
niteness, one may consider hard-core bosons to describe 
these bosons. Then, from such a bosonic Hubbard model 
one may derive the bosonic t-J model as its effective 
model for the case of strong on-site repulsion and small 
hole concentrations. This derivation is achieved just by 
following the steps developed in the theory of high-Tc su- 
perconductivity to derive the fermionic t-J model from 
the standard Hubbard model by tracing out the double- 
occupancy states. 

Therefore, to study the interplay of magnetic ordering 
and BEC of bosonic atoms, it is natural to start from 
the bosonic t-J model. In fact, very recently, Bonin- 
segni and Prokof'ev[l| studied this phenomenon by us- 
ing the bosonic t-J model where each bosonic electron 



is considered as a hard-core boson. They applied this 
model to study bosonic cold atoms 0, Hf. Because the 
model is purely bosonic, one can employ numerical anal- 
ysis. By quantum Monte Carlo(MC) simulations, they 
studied the low-temperature (T) phase diagram of the 
two-dimensional (2D) model for the case of anisotropic 
spin coupling Jx,y = aJz,a < 1, and found the coex- 
istence region of AF order and BEC as a result of the 
phase separation of hole-free and hole-rich phases. 

The bosonic t-J model has another important reason 
to study, i.e., it resembles to the fermionic t-J model of 
the high-temperature superconductors. There the inter- 
play of magnetism and superconductivity(SC) is of inter- 
est as one of the most interesting problems in strongly- 
correlated electron systems. At present, it is known that 
SC and antiferromagnetic (AF) Neel order can coexist in 
clean and uniform samples of the high-Tc cuprates [1]. 
Some theoretical works also report the coexistence of SC 
and AF order at T = o[5|. As this phenomenon appears 
as a result of interplay of fluctuations of quantum spins 
and BEC of superconducting pairs, simple mean-field-like 
approximations are inadequate to obtain reliable results 
to the relevant questions, e.g., whether the t-J model of 
electrons exhibits this coexisting phase. 

Not only these AF and SC transitions, the fermionic 
t-J model is expected to describe the metal-insulator 
transition(MIT) as observed in the cuprate high-Tc 
superconductors [I45I. At present, it seems that no well- 
accepted theoretical accounts for MIT beyond the mean 
field theory have appeared. 

In these situation, study of the bosonic t-J model be- 
yond the mean-field theory shall certainly shed some 
lights for understanding of these interesting problems in 
the fermionic t-J model. Then one can take advantage 
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of the bosonic nature of involved variables, which affords 
us to perform direct numerical simulations. 

In this paper, we study the bosonic t-J model. As 
explained above, our main motivations are the following 
two points; 

(a) studying the interplay of AF and BEC in the cold 
atoms; 

(b) getting insight for the AF, MIT and SC transitions 
of cuprate superconductors. 

We shall introduce a new representation of the bosonic 
t-J model for a s = ^ isotropic AF magnet with doped 
bosonic holes, and investigate the phase structure of the 
model by means of the MC simulations. Explicitly, we 
start with the slave-fermion representation of the original 
t-J model of electrons where the spinous are described by 
a CP^ (complex projective) field and holons are described 
by a one-component fermion field. Then we replace the 
fermion field by a Higgs field with fixed amplitude (U(l) 
phase variables) to obtain the bosonic t-J model. Thus 
the present model can be regarded as a slave-particle rep- 
resentation of the bosonic t-J model that is a canonical 
model for cold bosonic atoms in optical lattices. 

The usefulness of the slave-particle representation in 
various aspects has been pointed out for the original 
fermionic t-J model. We expect that similar advantage 
of the slave-particle picture holds also in the bosonic t- 
J model. One example is given by a recent paper by 
Kaul et al.Q. They argue a close resemblance in the 
phase structure between the fermionic t-J model and the 
bosonic one. Gapless fermions with finite density should 
destroy the Neel order. Furthermore, they induce a phe- 
nomenon similar to the Anderson-Higgs mechanism to 
the U(l) gauge dynamics, i.e., they suppress fluctuations 
of the gauge field strongly and the gauge dynamics is re- 
alized in a deconfinement phase [rj. Similar phenomenon 
is known to occur in the massless Schwinger model, i.e., 
(1-1-1) dimensional quantum electrodynamics, in which 
the long-range Coulomb interaction is screened by the 
gapless "electron" with a finite density of states Q. The 
Higgs field introduced in the present model to describe 
bosonic holons plays a role similar to these fermions, and 
then study of the bosonic t-J model may give important 
insight to the fermionic t-J model, in particular, proper- 
ties of low-energy excitations in each phase. 

The present paper is organized as follows. In Sec. 2, 
we shall introduce the bosonic t-J model in the CP^- 
spinon and Higgs-holon representation. We first consider 
the three-dimensional(3D) model at finite-T. An effec- 
tive action is obtained directly from the Hamiltonian of 
the fermionic t-J model by replacing the fermionic holon 
by bosonic Higgs holon. In Sec. 3, we exhibit the results 
of the numerical study of the model and phase diagram of 
the model. We calculated the specific heat, the spin and 
electron-pair correlation functions, and monopole den- 
sity. From these results, we conclude that there exists a 
coexisting phase of AF long-range order and the BEC in 
a region of low-T and intermediate hole concentration. 
In Sec. 4, we shall consider the two-dimensional (2D) sys- 



tem at vanishing T. We briefly review the derivation of 
the model. In Sec. 5 we exhibits various results of the 
numerical calculations and the phase diagram. Section 6 
is devoted for conclusion. 

II. BOSONIC t-J MODEL IN THE CP^-HIGGS 
REPRESENTATION: 3D MODEL AT FINITE T'S 

In this section, we shall introduce the bosonic t- 
J model in the CP^ spinon representation. In par- 
ticular, we first focus on its phase structure at finite 
temperature(T). To be explicit, let us start with the spa- 
tially three-dimensional (3D) original t-J model whose 
Hamiltonian H is given by 

H = -tY. {Ci+„.Cxa + H.c.) + • 1, 

Cxa — (1 ^ x(j^ '■'^(t) Cx(n 

Sx = ^C^aCx (ct : Pauh matrices), (2.1) 

where Cxa is the electron operator at the site x satisfying 
the fermionic anticommutation relations, a [ = 1(1), 2(J, 
)] is the spin index and 1 = 2,2 = 1 denote the opposite 
spin. /i(= 1,2,3) is the 3D direction index and also de- 
notes the unit vector. The first t term describes nearest- 
neighbor (NN) hopping of an electron without changing 
spin directions, i.e., ferromagnetic(FM) hopping. The 
second J term describes the nearest-neighbor(NN) AF 
spin coupling of electrons. The doubly occupied states 
{ClfCll\0)) are excluded from the physical states due 
to the strong Coulomb repulsion energy. The operator 
Cxcr respects this point. We adopt the slave-fermion rep- 
resentation of the electron operator Cxa as a composite 
form, 

Cxa = ^laxa, (2.2) 

where ipx represents annihilation operator of the 
fermionic holon carrying the charge e and no spin and a^o- 
represents annihilation operator of the bosonic spinon 
carrying s = 1/2 spin and no charge. Physical states 
satisfy the following constraint, 

CY o^xa^xa + V'lV'x)|phys) = |phys). (2.3) 

a 

In the salve-fermion representation, the Hamiltonian 
()2.ip is given as 

^ = ^l(^^x±i,0.xi>x±i, 

X + fJ. {a''aa)x. (2.4) 

We employ the path-integral expression of the parti- 
tion function Z = Trexp(— /3i/) of the t-J model in the 
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slave-fermion representation, and introduce the complex 
number a^g-^r) and the Grassmann number ipxi^) at each 
site X and imaginary time r G [0.J3 = {ksT)^^]. The 
constraint (|2.3p is solved [3 by introducing CP^ spinon 
variable Zxa-iT), i.e., two complex numbers Zxi,Zx2 for 
each site x satisfying 



and writing 



axa = (1 - 'ipx'ipxy^^Z., 



(2.5) 



(2.6) 



It is easily verified that the constraint (|2.3p is satisfied 
by Eqs. (|2.6p and (|2.5|) . Then, the partition function in 
the path-integral representation is given by an integral 
over the CP^ variables Zxai^) a-nd Grassmann numbers 

The bosonic t-J model in the slave-particle represen- 
tation is then defined at this stage by replacing ^px by 
a U(l) boson field (px (Higgs field) as 



ipx ^ \/S (px, (px = ex.p{i(px), 



(2.7) 



where S{— (ipl^px)) is the hole concentration per site 
(doping parameter), i.e., (^^Cl^Cxa) = 1 — S. Here 
we should mention that we have assumed uniform dis- 
tribution of holes and set the amplitude in front of (px 
a constant {\/S)- Validity of this assumption was partly 
supported by the numerical study of the closely related 
model with an isotropic AF coupling (ll. IT]|. Actually, it 
is reported there that the ground state of the bosonic t-J 
model (with Jx — Jy = Jz) is spatially uniform (without 
phase separation) for J/t < 1.5. We shall discuss on this 
point further in Sec. 6. 

We shall consider the system at finite and relatively 
high T's, such that the t— dependence of the variables 
Zx<n(px negligible (i.e., only their zero modes sur- 
vive). Then the kinetic terms of Zxa,(px including 
Zxdzx/dT,(pxd(px/ Bt disappear, and the T-dependence 
may appear only as an overall factor /3, which may be 
absorbed into the coefficients of the action and one may 
still deal with the 3D model instead of the 4D model. 

To obtain a description in terms of smooth spinon vari- 
ables, we change the CP^ variables Zxa at odd sites (the 
sites at which xi + X2 + xy, is odd) to the time-reversed 



CP variable 



,old 



as 



-old 



a' 

for ex = (-)^i+^2+^-^ = -1, 



(2.8) 



where e^fr' is the antisymmetric tensor e\i ~ —£21 = 
l,eii = £22 = 0. (Thus Zx\ = Zxi,Zxi = -Zx\) We 
stress that this is merely a change of variables in the 
path integral. (The AF spin configuration in the original 
variables becomes a FM spin configuration in the new 
variables.) 



In this way, the partition function Z of the 3D model 
at finite T's is given by the path integral, 

Z = J n [dzxd(PxY[dUx^.]expA, (2.9) 

x(xi,X2,X3) fJ- 

where the action A on the 3D lattice is given by 



A ^ As + Ai, 

Cl 

2 



— ~^ ^ ] (yZx+fJJxfiZx -\- C.C.^ , 



x^li 



^ ix^,-Zx4>x-,,(Px 
_xGeven,/j- 



^ ^ ^x + /i ^ 

a;Godd,/j, 



(2.10) 



where z' z = '^„z'^z„, etc.[l2[. We have introduced the 
U(l) gauge field Ux^j. = exp{i6x^) on the link {x,x -\- 
fi) {fi = 1,2,3) as an auxiliary field to make the ac- 
tion in a simpler form and the U(l) gauge invariance 
manifest. It corresponds to Ux^ ^ ZxZx+fj,/\zxZx+fj,\. 
In fact, one may integrate over Ux^ in Eq. (|2.9|) to ob- 
tain X^K/i fog^o(ci|^a:+/iZa;|) in the action (/q is the mod- 
ified Bessel function), which should be compared with 
the original expression J^xfj, \^x+^iZx\'^ ■ Both actions have 
similar behavior and it is verified for C3 = that they give 
rise to second-order transitions at similar values of ci . 

The action A is invariant under a local (x-dependent) 
U(l) gauge transformation. 



Zxa 
Uxfi 



e Zx 



(2.11) 



The gauge-invariant bosonic electron variable Bxa is ex- 
pressed as a composite as 



xa for ex = 1, 

xa for Ex = -1. 



(2.12) 



From Eqs. ()2.4p . ()2.10p . and the asymptotic form of 
log /q {x) , the parameters ci and C3 are related with those 
in the original t-J model as 



Cl 

C3 



J(3 for Cl >> 1, 

(2J/3)i/2 for Cl << 1, 

t5l3. 



(2.13) 



In terms of holons and spinons, one of our motivations 
is put more explicitly. In the slave-fcrmion representation 
of the fermionic t-J model, the holon-pair field ^xij. = 
ijjx'>Px+p, defined on the link behaves as a boson and a SC 
state appears as a result of its EEC. Hoppings of ^xf^ are 
to destroy the magnetic order because of their couplings 
to spinons as 



A 



SC 



x^fi^iy 



{Zx+i'Zx)^xti{Zx+tiZx+f_i+u)^x+u,fi, (2.14) 
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which is similar to Ah in Eq. (|2.10p . The difference is 
that (f>x is the site variable, while $3.^ is the link variable. 
Therefore, to study the CP^ model coupled with (fi^ may 
give us important insight to the dynamics of CP^ spinons 
and holon pairs on links, in particular the interplay of the 
AF order and SC. 

At the half filling {6 — 0), only the spinon part As 
survives, which describes the CP^ model. The parameter 
ci ~ J/{k-QT) = J 13 controls fluctuations of and Ux^i- 
The pure CP^ model As exhibits a phase transition at 
ci = ^ 2.8[ll|. In the low-T phase(ci > Cic), the 
0(3) spin variable, 

Sx = ZxOZx^ Sx ■ Sx — 1, (2.15) 

made out of Zxa has a long-range order, 
\\m\x\^oo{SxSa) — rn^ 7^ 0, i.e., the Neel order in 
the original model. (Note that the replacement (|2.8p 
leads to S'S™ = e^S"^'^, so an AF configuration of S°^'^ 
corresponds to a FM configuration of 5"°'". Throughout 
the paper we use the terms AF or FM configurations 
referring to the original spins S'^'^ .) The low-energy ex- 
citations are gapless spin waves, and the gauge dynamics 
is in the Higgs phase. The high-T phase(ci < Cic) is the 
paramagnetic phase, where the gauge dynamics is in the 
confinement phase and the lowest-energy excitations are 
spin-triplet bound states Sx of the spinon pairs. 

Let us next consider the role of the holon part A^. Ah 
comes from the hopping t term in (|2.1[) and the parameter 
C3 is expressed as in Eq. ()2.13l The spinon-pair amplitude 
of Zx and Zx+^ in As like Zx+^Zx measures the NN FM 
order of the original 0(3) Heisenberg spins due to the 
relation, 

\zx+^ • Zz\^ = -(1 — Sx ■ Sx+fi)- (2-16) 
Thus A\-i favors both a FM coherent hopping amplitude, 
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FIG. 1: Specific heat for L = 24 vs ci for fixed C3. It ex- 
hibits a peak of AF transition at c\ — 2.8(2.95, 3.25) for 
C3 = 0(0.5, 1.0). 
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FIG. 2: Size-dependence of the specific heat for C3 — 0.5. 

(zx+^Zx) and a BEC of (px- In other words, a BEG of 
4>x requires a short-range FM spin ordering. 



III. NUMERICAL RESULTS I (3D MODEL AT 
FINITE T'S) 

In this section we report the results of Monte Carlo 
simulations that we performed for the system given by 
Eas. (|2.9p and (|2.10p . We considered the cubic lattice 
with the periodic boundary condition with the system 
size (total number of sites) V = up to L = 30, and 
used the standard Metropolis algorithm. The typical 
statistics was 10^ MC steps per sample, and the averages 
and errors were estimated over ten samples. The typical 
acceptance ratio was about 50%. In the action we added 
a very small but finite (~ 10"'') external magnetic field. 

Let us start with the region of relatively high T's. In 
FiglHwe show the specific heat per site C , 

C^^{{A-{A)f), (3.1) 

as a function of ci for several values of C3. C exhibits 
a peak that shifts to larger ci as C3 is increased. This 
transition is nothing but the AF Neel phase transition 
observed previously for the C3 = case[l^. Fig [1] and 
the relations (|2.13p show that the Neel temperature is 
lowered by doped holes as we expected. 

In FigElwe present the system-size dependence of C for 
C3 = 0.5 in which the specific heat C develops moderately 
but systematically as we increase L. We fitted this C by 
using the finite-size scaling(FSS) of the form, 

C(ci,L) = L"/Xii/''e), (3.2) 

where t—{ci— cic)/cic with the critical coupling cic at 
infinite system size, and ct, v are critical exponents. In 
Figl3]we show the determined scaling function 0(x), from 
which we estimated the critical exponent of correlation 
length as = 1.7. 
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FIG. 3: Scaling function 
data of FigH 



(x) of p.2p determined from the 



In Fig [4] we present C as a function of C3 for fixed ci. 
There exist two peaks in C . To study the origin of these 
peaks, we define the "specific heat" Cs and Ch for each 
term As and Ah of the action (|2.10p separately as 



(3.3) 



In FigOwe present Cs and Ch for ci = 3.1, in which Cs 
has a peak at C3 ~ 0.8, the smaller-ca peak position of C 
while Ch has a peak at C3 ~ 1.3, the larger-ca position 
of C. From this result, we identify the peak at smaller 
C3 expresses the AF transition in Fig[Tl which is gener- 
ated by As, and the peak at larger C3 expresses the BEC 
transition, which is driven by Ah. Because each peak 
develops as L is increased, they are both second-order 
phase transitions. The critical exponent of Cs of Fig[5]is 
estimated as = 0.70. 

Let us turn to the low-T (large-ci) region and see what 



Cl=3.25 




0.5 _ 1.0 

Cs 



FIG. 4: Specific heat for L = 24 vs C3 for fixed ci. There are 
two peaks, which indicate the AF transition(in the smaller 
C3-region) and the BEC transition (in the larger C3 region). 
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FIG. 5: Specific heat Cg and Ch of (13.31) at ci = 3.1 and 
L = 24. The location of peak of Cs(h) almost coincides the 
location of peak of C in Fig|4]at smaller (larger) C3. 



happens to the AF and BEC phase transitions. In Figl6l 
we present C for ci — 6.0. Again there exist two peaks 
but the order of them is interchanged. Both peaks de- 
velop as L is increased, and therefore both of them are 
still second-oder phase transitions. As C3 (or equivalently 
i5) is increased, the transition into the BEC phase takes 
place first and then the AF phase transition follows. This 
means that there exists a phase in which both the AF and 
BEC long-range orders coexist. 

To verify the above conclusion, we measured the 
spin-spin correlation G^Fif) of Sx and the correlation 
Gbec {t) of the gauge-invariant "composite electron" pair 
variable V^^ in each phase. 



(3.4) 
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FIG. 6: C vs C3 for ci = 6.0 with L = 12, 18, 24, 30. There 
exist again two peaks of AF and BEC phase transitions. How- 
ever, their order are reversed (BEC peak in smaller C3 and AF 
peak at larger ca). Both peaks develop as L is increased, sup- 
porting that they are of second order. In the region between 
the two peaks, the AF order and the BEC order coexist. 
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FIG. 7: Spin-spin correlation function Gaf('") of p.4p in the 
three phases along ci — 6.0. For small C3, Gap(?") has a 
nonvanishing smooth component, i.e., the genuine AF stag- 
gered magnetization. (Remember the direction of the spin on 
the odd sites has been reversed.) For intermediate C3, there 
appears an oscillatory, i.e., FM component in the AF back- 
ground. For larger C3, SC exists without AF long-range order, 
where the FM long-range order dominates. 



GbEcW = (V^+r'^^f^Vxt,) +C.C., (3.5) 



v.. 



^x+fi^xj ^x — li 



^x+fj. 



^x ) ^x 1 • 



The results are shown in Figsl7]and[51 Fig [T] exhibits an 
interesting result in the coexisting phase of AF order and 
BEC at intermediate C3(oc 5). There the spin correlation 
has a FM component in the AF background. This shows 
that a short-range FM order is needed for the BEC as 
some mean-field theoretical studies of the t-J modelflj 
indicate. As C3 is increased further, the AF long-range 
order disappears and the FM order appears instead as a 
result of the holon-hopping amplitude. FiglS] shows that 
the BEC certainly develops for larger ca's. 

Let us comment here on the nature of BEC we have 
studied. To study BEC we have used Gbec, the 
long-range order of which implies a BEC of electron 
pairs. There may be another kind of BEC, a conden- 
sation of single bosonic electrons, which is measured by 
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FIG. 8: Correlation function Gbec('") (|3.5p of bosonic 
hole(vacant electron) pair 14^ in the three phases for ci — 6.0. 



FIG. 9: Monopole density p for ci = 6.0 and L = 24. Near 
the AF phase transition point at C3 = 4.8 (with an arrow), 
the curvature of p changes its sign. 



the electron-electron correlation function, {Bx+r,a-Bxa')- 
(Note that {(f)x+r4'x) is not suitable because it is not 
gauge invariant and always vanishes.) However, as we 
have seen in FigsmiH there is only one peak at most in 
the specific heat apart from the peak representing the 
magnetic order, so the two condensations, one of single 
electrons and the other of pairs of electrons take place 
simultaneously. 

We measured also the monopole density pfTs*!, which 
gives information on strongness of the fluctuation of the 
gauge field Uxf_i- In FigE] we present p along ci — 6.0. 
p changes its behavior around the AF transition point 
at C3 = 4.8, indicating that the AF ordered phase cor- 
responds to the deconfinement phase of Uxfj.- The low- 
energy excitations there are gapless spin waves described 
by the uncondensed component of Zx- In the paramag- 
netic phase, on the other hand, the confinement phase 
of Ux^ is realized and the low-energy excitations are the 
spin-triplet Sx with a gap. 



To summarize the results of the 3D system, we present 
in Fig llOl the phase diagram in the c^/ci — 1/ci, i.e., 

5 — T plane. The Neel temperature decreases slowly as 

6 increases, while the BEC critical temperature develops 
rather sharply. The two orders can coexist at low T's in 
intermediate region of 5. Fi gJTUl has a close resemblance 
to the phase diagram of Ref . However, in Ref. [l| , the 
anisotropic spin coupling is considered and the phase sep- 
aration occurs as a result. The coexistence phase of the 
AF and BEC orders obtained in Ref. [l| is nonuniform and 
accompanied with this phase separation. In the present 
paper, we studied the system with the isotropic spin cou- 
pling, and the phase with both the AF and BEC orders 
is realized under the uniform distribution of holes [15 1. 
In this sense, the model in the present paper is close to 
the fermionic t-J model with parameters t > J[5]. We 
expect that the fermionic t-J model has a similar phase 
diagram to Fig |lojT6j . In Sec. 6, we explain further the 
implication of the results obtained in the present paper 
to the phase structure of the fermionic t-J model. 
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FIG. 10: Phase diagram of the 3D model at finite T in the 
C3 — ciplane (a) and (b) in the cs/ci — 1/ci plane (cs/ci oc tS/J 
and 1/ci oc T). In the low-T and intermediate 5 region, there 
exists the coexisting phase of BEG and AF order. All the 
phase transition lines are of second order. 



IV. BOSONIC t-J MODEL IN THE CP^-HIGGS 
REPRESENTATION: 2D MODEL AT T = 



As mentioned in Sect.l, the cold atoms can be put on 
a 2D optical lattice. At T = one may expect a BEC 
and/or magnetic ordering at certain conditions for den- 
sity of atoms per well, interaction between atoms, etc. 
We expect that the bosonic t-J model in the present sec- 
tion describes dynamics of bosonic atoms with (pseudo- 
)spin degrees of freedom and (strong) repulsive interac- 
tions between them. For the physics of high-Tc cuprates, 
it is also interesting to study the two-dimensional {2D) 
bosonic t-J model at T = 0. The main concern there 
is the quantum phase transitions (QPT) of magnetism, 
MIT, and SC. Study of the QPT in the doped CP^ model 
is not only important for verifying the phase diagram at 



finite T obtained in the above, but also interesting for 
physics of cold atom systems in optical lattices. 

In this section, we shall briefly review the derivation 
of the effective field-theory model for the 2D fermionic 
t-J model at T = Ofg?]. Then, in this effective model, by 
replacing the fermionic holon variables by Higgs boson 
variables, we obtain the bosonic effective model of the 2D 
bosonic t-J model at T = 0. The main difference from the 
previous sections (the 3D model at finite T's) is that the 
action contains the kinetic terms of spinous and holons 
due to their nontrivial T-dependence, and one should con- 
sider the 2+1 dimensional model. They are both 3D 
models but the couplings along the third-direction (2:3 or 
r) are different in the two models. 

Let us start with the path-integral representation of 
the (quantum) partition function Zp^ of the 2D fermionic 
t-J model in the slave-fermion representation with the 
CPi variables 0, 



r,2D 



J l[[dMr)dzr{T)] exp (a|D(t)) , (4.1) 



where r = {xi,X2) denotes the site of the 2D lattice and 
r G (0, 00) is the imaginary time. The action A'^{t) is 
given bv[l3| 



f°° r 

/ dr'S^ — [iprDrr'^r + ZrdrZr) 

Jo ^ L 



DrT — 



1=1,2 



4 ^ ' 

i 



(4.2) 
(4.3) 



where i ~ 1,2 denotes the spatial direction index and 
also the unit vector. 

Each kinetic term, ■iprDrT'4'r, ZrdrZr, in is purely 
imaginary, so the straightforward MC simulations cannot 
be applied. However, concerning to ZrdrZr, by integrat- 
ing out the half of the CP^ variables (e.g., z^'s on all 
the 2D odd sites, i.e., ef^ = (-)^i+^2 ^ _i) one may 
obtain a purely real action as a result. That real action 
describes the low-energy spin excitations in a natural and 
straightforward manner. 

Explicitly, we assume a short-range AF order, 



{Z(7z)r ^ —(zaz) 



(4.4) 



Then we parameterize each CP^ variable Zq at the 2D 
odd site o by referring to one of its NN partner z^ at the 
even site e at equal time as 



Zq — PocZq + (1 IT^OC I ) ^ UqqZqj 



(4.5) 



where Poe is a complex number sitting on the link (oe) 
and Uoe is a U(l) phase factor that makes the param- 
eterization (14. 5p to be consistent with the local gauge 
symmetry [17|- The assumption (j4.4p implies |pocP ^ 1, 
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which is justified by the J-term in the action (|4.2p for 
the hghtly doped case {S <C Integration over Zq is 

reduced to the integral over poe, which may be approxi- 
mated as the Gaussian integration Unk by Hnk, 



dpdp exp(— pMp + pk + £p) 
oc detM"^ exp(£M-^fc). (4.6) 



The CP^ part of the resulting action is not symmetric 
concerning to spatial directions because the choice of a 
definite even-site partner breaks it. However, by consid- 
ering the smooth configurations of spinous 2;(r, r), one 
recovers the symmetric action in the form of continuum 
space (T^, 



A, 



2a2 



dr / d^r 



X 



Ja^ DizDiZ + —DT-zDrz\ , 



2J 

1=1,2 

= iA^, Af_, = izdf_,z, 



(4.7) 



where a is the lattice spacing (hereafter we often set a = 
1). The relativistic couphngs of (|4.7p give rise to spin 
wave excitations in the background of AF order with the 
dispersion E{k) oc k instead of E{k) oc fc^ in the FM spin 
model. 

The integration over the odd-site spins affects the 
holon and spinon hopping term. The most important 
point is that the U(l) factor Uoc always appears in com- 
bination with Tpo as Uoe^^o, and therefore we redefine 
C^ooV'o V'o- Then new ipo transforms under a gauge 



transformation as "00 



tpo- This property is con- 



sistent with the fact that spinon hopping amplitude, e.g., 
from even site (e') to odd site (o), contains factor [z^-z^:) 
in the resultant effective action like 



(Ze • Zc-)lpotpc' + C.C. 



(4.8) 



Finally, the holon kinetic term is now expressed by 
using Vo™ = C/ocC''* aslU 



dr 



E 



(4.9) 



To proceed we note that the obtained total action is not 
symmetric w.r.t. the spatial directions, but a symmetric 
one can be available by reintroducing CP^ variables at 
the 2D odd sites according to the invariance under a lo- 
cal U(l) gauge transformation. Furthermore, by choosing 
the lattice spacing in the direction suitably, the sym- 
metry in all the three directions is recovered. Then the 
spin part of the action becomes nothing but the 3D 
action of (|2.10p . Both of them have the same con- 
tinuum limit (j4.7p for the spin part, so they are to be 
categorized to the same universality class. 



Then the bosonic t-J model on the lattice is obtained 
by the replacement ipx \/54>x^ 4>x = s^piiy^x) as 
in Eq. ([2Jl) . The partition function Z^^ on the 2+1- 
dimensional spacetime lattice is then given by 



Z 



2D 



A'"" - 

— 



Y[[dz,d(b^ Y[ dU^^] exp(A2D), 
+ A^D, = As of Eq.dllOl), 

Y{^^^+3Ux3't'x + ^4>x+3Ux34'x (4-10) 

x^c x^o 

Zx+iZx(t>x+i(t>x + ZxZx+i4'x+i4>x + C.C.^ . 



x^G,i 



Here x = (a;i, a;2, xa), and r = (xi,a;2) stands for the 
coordinates of the 2D plane and X3{= 0, 1, • • • , oo) is the 
discretized imaginary time (r — x^a^). Even (e) and 
odd (o) sites are defined regarding to the 2D plane as 
g2D ^ ^_-^x,+x2 = i(even), -l(odd). 

The coefficients ci , S3 are related to the parameters of 
the t-J model as follows: 1231 



ci constant independent of J and <, 

S3 - jS. 



(4.11) 



ci measures the solidity of the Neel state in the AF 
Heisenberg magnet at T = 0. S3 is the hopping am- 
plitude of electrons. For the AF Heisenberg model with 
the standard NN exchange coupling, ci > cic — 2.8, and 
the critical value cic decreases (increases) if long-range 
and/or anisotropic couplings that enhance (hinder) the 
AF order are added [2^. 

In All we added the Hermitian conjugate of the holon 
kinetic term, (t>x+3(Ux3,Ux3)4'x, which were absent in 
Ea. (|4.9p . to make the action real, explicitly. This mod- 
ification might sound crucial. However, we expect it a 
modest replacement because the omitted imaginary part, 
~i{(j)x+3Ux3(f>x—C-C.), etc, should have vanishing expecta- 
tion value and behave mildly for the relevant configura- 
tions for the original action without the complex conju- 
gates. Another support for this modification is given by 
taking into account the fiuctuation of holon field [24j. 

Let us summarize the difference between the 2D model 
at T = and the 3D model at T > of ([TTU)) . To see it 
explicitly, it is convenient to use the redefinition. 



{ 



for 
for 



-1 



(2D model), (4.12) 



and the relation zz' ~ —zz'. Then A^ becomes 

X 



In the similar notation, 

, , ^ r 4>x for ex 
<t>x for ex 



(3D model). 



(4.14) 
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the action Ah of the 3D model (|2.10p is given as 

= y ^ (ea;Z2:+MZa;(/>U^</>^ + C.C.y (4.15) 

Thus, the hopping of holons (l)^ along fi — 3 is uniform 
in the (/i0-channel in the 2D model, while it is alternative 
in the 0(/)-channel in the 3D model, i.e., the latter has an 
extra factor (— )^^. 



V. NUMERICAL RESULTS II (2D MODEL AT 

r = o) 

A. Phase diagram 

We studied the quantum system A^^ of (|4.10p by MC 
simulations ^5] . In this section we show the results of 
these calculations. Let us first present the obtained phase 
diagram in Fig llll Its phase structure is globally simi- 
lar to Fig[TO] for the 3D model at finite T's. As S3 is 
increased, the AF order is destroyed and the BEC ap- 
pears. The appearance of a first-order transition line is 
new. Below we shall see the details of various quantities, 
which support Fig[Tl] We select the following four points 
in the phase diagram to represent each of four phases; 

(I) S3 = 1, ci = 2; Normal and paramagnetic state. 

(II) S3 = 1, ci = 20; AF state without BEC. 

(III) S3 6, ci = 4; BEC state with a FM order 

(IV) S3 = 30, ci = 20; AF and BEC state. 

As presented later, various quantities are measured on 
these points and compared with each other. 

Let us comment here on the assumption of the short- 
range AF order (|4.4p in deriving the effective model 
(|4.10p . As the phase diagram of FigfTTl shows, we find 
that there exists a FM order for small ci and large S3. So 
the applicability of (j4.4p in this parameter region might 
be questionable. However, the existence of the FM or- 
der itself in Fig[Tl] is consistent with the results of the 
3D finite-T system obtained in Sees. 2, 3 (as long as each 
phase at T = smoothly continues to T > 0). Thus we 
expect that the global phase structure of FigHT] is not 
affected by the assumption (|4.4p . although the location 
and the nature of the phase boundaries need some mod- 
ifications. Furthermore, calculations of various physical 
quantities in that phase give us very important insight in 
understanding the whole phase structure of the model. 



B. Specific heat and internal energy 

We first study the "specific heat" (fiuctuation of the 
action A^^), 

C2D = 1((A2D_(a2D))2). (5.1) 



AF 




0.0 0.5 1.0 1.5 2.0 
S3/C1 

(b) 



FIG. 11: Phase diagram of the 2D t-J model at T = in the 
S3 — ci plane (a) and (b) in the sa/ci — 1/ci plane. The ci-term 
controls the stability of the AF Neel order, whereas the S3- 
term represents holon hopping. In the large-ci and interme- 
diate S3 /ci , there exists the coexisting phase of AF order and 
BEC. At the merging point of two transition lines, the peak 
in C^^ disappears. NS stands for normal state and FM for 
ferromagnetic order. The phase transition line separating AF 
and BEC-I-AF phases is of first order, while other three tran- 
sition lines are of second order. The four points (I-IV) marked 
by stars are selected as typical points: (I) ci — 2,8^ — 1; (II) 
ci = 20, S3 = 1; (III) ci = 4, S3 = 6; (IV) ci = 20, S3 = 30. 



In Fig[T2]we present C^^ as a function of S3 for various 
values of ci. We also present the "specific heats" C^^ 
and Cf^ of each term Af-* , Af^ in the action as defined 
in (|5.ip . As in the previous 3D finite-T case, these results 
indicate that there exist two phase transition lines in the 
S3 — C1 plane. They intersect at ci ~ 5.0, S3 ~ 4.0. At the 
BEC transition points at high Ci 's (low 1/ci ) (see FigfT^ 
for ci=8.0 near S3 ~7.0), C^^ exhibits large values and 
large system-size dependence. 
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FIG. 13: "Internal energies" (A^^), (A^^), {Ai°) vs 1/ci for 
fixed S3/C1 = 0.83. We also show C^^ . The intersection point 
of the two phase transition lines (merging point of the AF and 
BEG phase transitions) is marked by an oval. At this region, 
C^^ exhibits no anomalous behavior. Slightly anomalous be- 
havior of (Ag^) and (A^) cancel with each other and the 
total {A^^) is a regular function of 1/ci. 




6.5 7.0 7.5 S3 8.0 8.5 9.0 



FIG. 12; "Specific heats", C^" , C^^" , vs S3 for ci = 
4.0, 5.0, 8.0. The notation Cs - 24 denotes for L = 24, 
etc. At Cl = 5.0 the two peaks almost merge. 



In order to see what happens at this intersecting point 
of the two phase transition lines, ci ~ 5.0, S3 ~ 4.0, we 
measured "internal energies" (A^^), {A'^^), (A^), and 
C^D, C2D, as a function of 1/ci along 33/01 = 0.83. 
We show the result in Fig[T31 At the intersection point, 
the total specific heat C^^ exhibits no anomalous be- 
havior, though Cg^ and show peaks at that point. 
The calculations shown in Fig |13l show an interesting phe- 
nomenon that anomalous behavior in {A'^^) and (A^) 
cancel with each other and the total {A"^^) is a regular 
function of 1/ci. 

To verify the order of the phase transitions, we mea- 
sured distribution of A"^^ for configurations generated 
through the MC steps. In Fig[T3]we present the distribu- 
tions of around ci = 12, S3 = 10, 11, 12. The distri- 
bution at the middle point, S3 — 11, shows a double-peak 
structure, so we concluded that the phase transition at 
this point is of first order. We found that the phase tran- 
sition separating the AF and AF-I-BEC phases is of first 
order, while the other three transition lines are of second 
order. The difference of the order (1st vs. 2nd) in the 
3D and 2D models may be attributed to the difference of 
holon hopping along the 3rd direction explained at the 



end of Sect. 4. The present 2D model has a (pcf) uniform 
coupling, which put more weight for ordering of holons 
and spinons than the 3D model. This may generate the 
first-order transition in certain regions. 




- 0^ = 12,53=12 




FIG. 14: Distribution of A^° . At ci = 12, S3 = 11, it exhibits 
a double-peak structure, whereas the other two do not. This 
result means that the phase transition here is of first order. 
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C. Spin correlations 

In this subsection, we investigate the spin correlation 
function at equal time, 

GafW= (5.2) 

and show the snapshots of spin configurations. In Fig [15] 
we present G'af('') and the snapshots of spins at the 
typical four points (I-IV) in the phase diagrams of 
FigHH In the snapshots, each Sx starts from the center 
of the unit sphere and ends at a dot on the sphere. 

(I) ci = 2, S3 = 1 : Directions of spins are random, and 
there is no long-range magnetic order. 

(II) ci = 20, S3 = 1 : This point in the phase diagram is 
located in the deep AF region. It is obvious that there 
exists the AF long range order. 

(III) ci = 4, S3 = 6 : There appears the oscillative order 
around Gaf — 0, i.e., the FM long range order instead 
of the AF order. 

(IV) ci = 20, S3 = 30 : The even and odd site spins 
have their own magnetizations, and Mq, and these 
even and odd magnetizations cant with each other. This 
corresponds to the canting state studied in Ref.[14]. In 
other word, there appear a component of a FM order in 
the background of AF long-range order. 




-1.0-U , — . — , — . — , . — ,-r -1.0-U — , — . — , — . — , — . — ,-r 

4 8 12 16 4 8 12 16 



FIG. 15: Spin correlation function Gaf('") of (|5.2|) and spin 
snapshot at the points (I-IV) in Fig llll (I) ci =2, S3 = 1. 
There is no long-range order of spins. (II) ci — 20, S3 = 
1. Endpoints of spins are centered near the south pole of 
the sphere, indicating the (AF) long-range order. (Ill) ci = 
4, S3 = 6. Results indicate the alternative order, i.e., the FM 
order. (IV) Ci = 20, S3 = 30. Endpoints form two groups on 
the sphere, indicating a canting "order" of spins [14j. 



D. Gauge dynamics and topological objects 

In this subsection, we study the gauge dynamics asso- 
ciated with two U(l) gauge fields; the spin hopping am- 
plitude Uxfi and the dual U(l) holon hopping amplitude 
Wxi defined as follows; 

^ e U{1). (5.3) 

For Uxn^ we calculated the instanton density p as in 
Sec. 3 [26], whereas for Wxi, we calculated its gauge- 
invariant flux density, 

= 2^ (mod [~i\og{Wx2Wx+2.iWx+i,2Wxi) ,2Tr]) , 

(5.4) 

as Wxi has only the spatial components. In Fig[Tn] we 
present p and fw as functions of S3; (a) ci = 8.0 and (b) 
Ci = 4.0. For the both cases, it is obvious that p and 
fw change their values suddenly at the relevant phase 
transition points determined by C^^. We also verified 
that the singular behavior of the "specific heat" Cg (Ch) 
correlates to that of p (fw) as it is expected. 

p influences the low-energy excitations of the spinon 
sector, whereas fw is related to the holon hopping. The 
large-p region corresponds to the confinement phase of 
spinons and the low-energy excitations there are de- 
scribed by the spin-triplet Zx(JZx. On the other hand, 
in the small-p region of the AF state, deconfinement of 




FIG. 16: The instanton density p and the flux density fw 
of ((El for (a) ci = 8.0 and (b) ci = 4.0. "TP,"("TPh") 
indicates the location of the AF (BEG) phase transition point. 
p changes suddenly at TPs. while fw changes suddenly at 
TPh as expected. 
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quanta takes place, and the low-energy excitations are 
gapless spin waves. 

Large fluctuations of W^i hinder coherent hopping of 
holons, and induce large fluctuations of the holon field (\)x- 
However, as the holon density 8 is increased, the holon 
hopping term stabilizes Wxi and reduces jw ■ From our 
study of p and /vk, there seem to be some correlations 
between them, the origin of which is of course the term 
AT. 



E. Electron correlations 

In this subsection, we study the correlation functions 
of single electrons and electron pairs in the 2D lattice 
in the various phases. We define the bosonic electron 
operator, B^a^ as follows; 



Bxa is invariant under the local U(l) gauge transforma- 
tion. We define also operators of the spin-singlet and the 
spin-triplet pairs of bosonic electrons on the NN sites at 
equal time. 



AS _ Ht Rt _ kA K?t 

A* — + 



(5.6) 



Then we introduce correlation functions of B^a and A^'^ 
at equal time as 

x.i.a 

^''*(^) ^ ^ E ^Xr^,^x^■ (5.7) 
x,i^j 

Their behaviors in the previous four points (I-IV) in 
FigHDare shown in Fig[T71 Each point has the following 
properties; 

(I) NS phase and (II) AF phase: All the three functions, 
GB{T)iG'^"^iy) have no long-range correlations. 

(III) BEC-fFM phase: G'B(r) and the triplet pair G'*(r) 
exhibit the long-range order, and the BEC takes place in 
these channels. 

(IV) BEC-I-AF phase: all the three functions exhibit the 
long-range order, and the BEC takes place. 

As discussed in Sect. 3, these results are consistent with 
the expectation that the two BECs, one of single elec- 
trons and the other of (triplet) electron pairs, take place 
at the same time. Furthermore, these results and the 
previous results on the (AF and FM) spin orders indi- 
cate that the BEG order and spin order can be superim- 
posed in certain region. This is another example of the 
"spin-charge" separation. 
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FIG. 17: Electron correlation functions of (|5.7p at the four 
points (I-IV) of Fig fm (a) G'fl(r) of single bosonic electrons; 
(b) G°''(r) of spin-singlet and triplet pairs. Their meaning in 
each point is discussed in the text. 



VI. CONCLUSION 

In this paper, we have investigated the bosonic t-J 
model in the CP^-spinon -f U(l) Higgs-holon representa- 
tion. This model is introduced by replacing the fermionic 
holons with the bosonic ones in the slave-fermion t-J 
model. We studied its phase structure and the dynami- 
cal properties both in the 3D finite-T model and the 2D 
T — model. In particular, we are interested in the in- 
terplay of the BEC and AF order, because a coherent 
holon hopping is required for the BEC, whereas it hin- 
ders the AF long-range order. Our study by means of 
the MC simulations exhibits a phase diagram in which 
the coexistence region of AF and BEC orders appears. 
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This result suggests that in the fermionic t-J model a 
MIT takes place at a finite hole concentration 5^ and 
that phase transition point is located within the AF re- 
gion of the spin dynamics. From the result of the present 
paper, we also expect that as 5 is increased further, the 
BEC phase transition through the formation of electron 
pairs and their BEC takes place in the AF region. Ac- 
tually, this problem can be studied in the framework of 
the fermionic t-J model in the slave-fermion representa- 
tion. Explicitly, one may obtain the effective model by 
integrating out the fermionic holon field by the hopping 
expansion. The resultant model includes only the bosonic 
variables; Zx of spinous, "order parameter" of the coher- 
ent hopping of holons (i.e., the dual gauge field Wxi)^ and 
the SC order parameter of electron pairs, which can be 
analyzed numerically (1^. For this scenario one should 
include some interaction terms that we have ignored in 
Sect. 4 in our action, but have been generated in the pro- 
cess of integrating half of the spinon variables [9)]. Among 
these terms, there is an attractive interaction between 
fermionic holons at NN sites [l^l- This term favors for- 
mation of holon pairs, and their coherent condensation 
gives rise to the superconductivity. 

In an optical lattice, cold atoms at each well can have 
(pseudo-)spin degrees of freedom 2s -I- 1 = 1, 2, • ■ • , and 
there appears an AF interaction among them as a re- 
sult of the repulsion. In this case, the CP^ constraint is 
changed to X]a=i l-^^xal^ ~ 2s, but this change of the nor- 
malization of Zx can be absorbed by the parameter ci. 



Then cold atom systems with one particle per well can 
be described by the CP^ model hke of Eg . (prU)) p8| . 

We expect that A^^ of (|4.10p describes a 2D system of 
bosonic cold atoms with (pseudo-)spin and repulsive in- 
teraction in an optical lattice slightly away from the case 
of one-particle per well. For such system of cold atoms 
the phase diagram, Fig[Tll indicates the coexistence of a 
BEC and a long-range order of an internal symmetry 0] . 

Let us comment on the effect of our replacement (|2.7p 
of fermionic holons by Higgs field with a uniform ampli- 
tude, y^(/)x- Being compared with the faithful bosonic 
model of the original t-J model, this replacement cer- 
tainly ignores fluctuations of the amplitude (density) of 
holon variable, which disfavors the BEC of holons. How- 
ever, we expect that the BEC we obtained in the present 
model should survive in the faithful bosonic t-J model, 
although location of the transition curves may change. 
This problem is under study and we hope that result will 
be reported in a future publication [2^. 
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